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We analyze the dynamics of entanglement in a two-qubit system interacting with its environment 
via a quantum nondemolition system-reservoir interaction. The system and reservoir are initially 
assumed to be separable with the reservoir being in an initial squeezed thermal state. Since the 
resulting dynamics becomes mixed, in order to study the ensuing entanglement, we make use of a 
recently introduced measure of mixed state entanglement via a probability density function which 
gives a statistical and geometrical characterization of entanglement by exploring the entanglement 
content in the various subspaces spanning the two-qubit Hilbert space. After developing the general 
dynamics, which is categorized into the independent or collective decoherence regime, of N qubits 
interacting with their bath (reservoir) via a quantum nondemolition interaction we specialize to 
the two-qubit case for applications. The dynamics is found to satisfy hermiticity and spin-flip 
symmetry operations. We give an effective temperature dependent dynamics for the collective 
decoherence regime and make an application of the two-qubit dynamics to a simplified model of 
quantum repeaters. 

PACS numbers: 03.65.Yz, 03.67.Mn, 03.67.Bg, 03.67.Hk 



I. INTRODUCTION 



Open quantum systems are ubiquitous in the sense that any system can be thought of as being surrounded by its 
environment (reservoir or bath) which influences its dynamics. They provide a natural route for discussing damping 
and dephasing. One of the first testing grounds for open system ideas was in quantum optics Its application 
to other areas gained momentum from the works of Caldeira and Leggett 0, and Zurek among others. The 
total Hamiltonian is H = H$ + Hr + Hsr , where S stands for the system, R for the reservoir and SR for the 
system-reservoir interaction. The evolution of the system of interest S is studied taking into account the effect of 
its environment R, through the SR interaction term, making the resulting dynamics non-unitary. Depending upon 
the system-reservoir (S* — R) interaction, open systems can be broadly classified into two categories, viz., quantum 
non-demolition (QND) or dissipative. A particular type of quantum nondemolition (QND) S — R interaction is given 
by a class of energy-preserving measurements in which dephasing occurs without damping the system, i.e., where 
[Hs,Hsr] = while the dissipative systems correspond to the case where [Hs,Hsr] ^ resulting in decoherence 
along with dissipation Q . 

A class of obscrvables that may be measured repeatedly with arbitrary precision, with the influence of the mea- 
surement apparatus on the system being confined strictly to the conjugate observables, is called QND or back-action 
evasive observables [H, B, B, B ■ Such a measurement scheme was originally introduced in the context of the detection 
of gravitational waves [3, [lfj. The energy preserving measurements, referred to above, form an important class of 
such a general QND measurement scheme. 

The interest in the relevance of open system ideas to quantum information has increased in recent times because 
of the impressive progress made on the experimental front in the manipulation of quantum states of matter towards 
quantum information processing and quantum communication. Myatt et al. and Turchette et al. fl2| have 

performed a series of experiments in which they induced decoherence and decay by coupling the atom (their system-,!?) 
to engineered reservoirs, in which the coupling to, and the state of, the environment are controllable. An experiment 
reported in Ref. (]~3j | demonstrated and completely characterized a QND scheme for making a nondeterministic 
measurement of a single photon nondestructively using only linear optics and photo-detection of ancillary modes, to 
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induce a strong nonlinearity at the sing le photon level. The dynamics of decoherence in continuous atom-optical QND 
measurements has been studied in [14| . 

Quantum entanglement is the inherent property of a system to exhibit correlations, the physical basis being the 
non-local nature of quantum mechanics [15j . and hence is a property that is exclusively quantum in nature. Entan- 
glement plays a central role in quantum information theory [l6| as in interesting non-classical applications such as 
quantum computation |17l | and quantum error correction (l8j | . A number of methods have been proposed for creating 
entanglement involving trapped atoms [l^, Hl[ . 

An important issue is to study how quantum entanglement is affected by noise, which can be thought of as a 
manifestation of an open system effect |22| . A recent experimental investigation of the dynamics of entanglement 
with a continuous monitoring of the environment, i.e., via a realization of quantum trajectories p3j . has been made 
in (24|. Here we study the effect of noise on the entanglement generated between two spatially separated qubits, 
by means of their interaction with the bath, which is taken to be in an initial squeezed-thermal state [J, [25[ . This 
is of relevance to evaluate the performance of two-qubit gates in practical quantum information processing systems. 
The two qubits are intially uncorrelated. With the advent of time entanglement builds up between them via their 
interaction with the bath but eventually gets destroyed because of the quantum to classical transition mediated by 
the noise. In this paper we study the effect of noise generated by a QND S — R interaction and take up the issue of 
a dissipativc noise in a separate work. 

Since we are dealing here with a two qubit system which very rapidly evolves into a mixed state, a study of 
entanglement would necessarily involve a measure of entanglement for mixed states. Entanglement of a bipartite 
system in a pure state is unambigious and well defined. However, mixed state entanglement (MSE) is not so well 
defined. Thus, although a number of criteria such as entanglement of formation [26|, [2?], [28[ and separability [2!| exist, 
there is a realization [261 ] that a single quantity is inadequate to describe MSE. This was the principal motivation for 
the development of a new prescription of MSE [30j | in which it is characterized not as a function, but as a probability 
density function (PDF). The known prescriptions such as concurrence and negativity emerge as particular parameters 
that characterize the probability density. We will principally make use of this measure in our study of entanglement 
in the two-qubit system. 

The plan of the paper is as follows. In Section II, we give a brief description of the recently developed entanglement 
measure of MSE [301 ] . In Sections III and IV, we develop our open system model for the multi-qubit dynamics under 
the influence of a QND S — R interaction. Section III develops the general dynamics for a multi-qubit system, where 
the qubits are spatially separated, initially uncorrelated and the bath is in a general squeezed-thermal state. Section 
IV specializes these considerations to the case of two qubits. In Section V, we point out some interesting symmetries 
obeyed by the two-qubit dynamics. Section VI deals with the entanglement analysis of the two-qubit open system 
using the PDF as a measure of entanglement. We also briefly dwell upon the usual measure of MSE, concurrence. We 
deal with the scenarios where the two qubits effectively interact via localized S — R interactions, called the independent 
decoherence model, as also when they interact collectively with the bath, called the collective decoherence model. In 
Section VII, we make a brief discussion of the temperature dependent effective dynamics obeyed by the two-qubit 
open system in the collective decoherence regime. Section VIII makes a brief application of the model to practical 
quantum communication, in particular, in the realization of a quantum repeater |3lll32T]. In Section IX, we make our 
conclusions. 



II. CHARACTERIZATION OF MIXED STATE ENTANGLEMENT THROUGH A PROBABILITY 

DENSITY FUNCTION 

Here we briefly recapitulate the characterization of mixed state entanglement (MSE) through a PDF as developed 
in [30( | . As pointed out in the Introduction, the above criterion was evolved from the motivation that for the charac- 
terization of MSE, a single parameter is inadequate. The basic idea is to express the PDF of entanglement of a given 
system density matrix (in this case, a two-qubit) in terms of a weighted sum over the PDF's of projection operators 
spanning the full Hilbert space of the system density matrix. The PDF of a system in a state which is a projection 
operator p = jj^-ai of rank M is defined as: 

rn M {t) = T~ru ' W 

J aH-iiM 

where f dH]j M is the volume measure for Hn M , which is the subspace spanned by IIm. The volume measure is 
determined by the invariant Haar measure associated with the group of automorphisms of J dTi]j M , modulo the 
stabilizer group of the reference state generating 7iu M ■ Thus for a one dimensional projection operator, representing 
a pure state, the group of automorphisms consists of only the identity element and the PDF is simply given by the 
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Dirac delta. Indeed, if p — III = the PDF has the form V p {£) = 8(£ — £^,) thereby resulting in the description 

of pure state entanglement, as expected, by a single number. The entanglement density of a system in a general mixed 
state p is given by resolving it in terms of nested projection operators with appropriate weights as 

p = (Ai - A 2 )n a + (A 2 - A 3 )n 2 + (Ajv_i - A JV )n A r„ 1 + X N U N 

N 

= A M^ M , (2) 

M=l 

where the projections are U M = YljLi IV'jXV'j'l) w ^ tn M = 1,2, ...,7V and the eigenvalues Ai < A 2 < i.e., the 
eigenvalues are arranged in a non-decreasing fashion. Thus the PDF for the entanglement of p is given by 

N 

V p (£) = u M Vn M (£), (3) 

Af=l 

where the weights of the respective projections Vn M (£) are given by lum = Ajvf/Ai. For a two qubit system, the density 
matrix would be represented as a nested sum over four projection operators, IT, II2, IT3, II4 corresponding to one, 
two, three and four dimensional projections, respectively, with 111 corresponding to a pure state and II4 corresponding 
to a a uniformly mixed state, is a multiple of the identity operator. The most interesting structure is present in II2, 
the two-dimensional projection, which is characterized by three parameters, viz. £ C usp, the entanglement at which 
the PDF diverges, £ max , the maximum entanglement allowed and V<x{£ ma ^), the PDF corresponding to £ max . The 
three dimensional projection II3 is characterized by the parameter £±, which parametrizes a discontunity in the 
entanglement density function curve. By virtue of the convexity of the sum over the nested projections ([2]), it can be 
seen that the concurrence of any state p is given by the inequality C p < (Xi — A 2 )Crii + (^2 — A 3 )Cri2- Thus while the 
concurrence for a three and four dimensional projection is identically zero, through the PDF one is able to make a 
statement about the entanglement content of these spaces. Also, as pointed out in [3(|, in the case of NMR quantum 
computation, concurrence and negativity are zero, whereas the PDF is able to elucidate the role of entanglement 
utilized by the NMR operations. These features as well as the fact that the PDF (J3j) enables us to study entanglement 
of a physical state by exploiting the richness inherent in the subspaces spanned by the system Hilbcrt space makes 
the PDF an attractive statistical and geometric characterization of entanglement. We provide an explicit illustration 
of this in Section VI. 



III. TWO-QUBIT QND INTERACTION WITH A SQUEEZED THERMAL BATH 

We consider the Hamiltonian, describing the QND interaction of L qubits with the bath as 0, [H, [34| 

H = Hs + Hr + Hsr 

L 

= £ Ke n J: + Y tokblh + E hJ z(9kb{ + gTh). (4) 

n k n.k 

Here Hs, Hr and Hsr stand for the Hamiltonians of the system, reservoir and system- reservoir interaction, respec- 
tively, fet, bk denote the creation and annihilation operators for the reservoir oscillator of frequency u>k, 5^ stands 
for the coupling constant (assumed to be position dependent) for the interaction of the oscillator field with the qubit 
system and are taken to be 

9t = 9ke- ik - r % (5) 

where r n is the qubit position. Since [Hs, Hsr] = 0, the Hamiltonian (1) is of QND type. In the parlance of quantum 
information theory, the noise generated is called the phase damping noise [1, [35[ . 

The position dependence of the coupling of the qubits to the bath (f5|) helps to bring out the effect of entanglement 
between qubits through the qubit separation: r mn = r m — r n . This allows for a discussion of the dynamics in two 
regimes: (A), independent decoherence where k.r mn ~ ^j^- > 1 and (B). collective decoherence where k.r mn ~ ^j^- — > 
0. The case (B) of collective decoherence would arise when the qubits are close enough for them to experience the 
same environment or when the bath has a long correlation length (set by the effective wavelength A) compared to the 
interqubit separation r mn Our aim is to study the reduced dynamics of the qubit system. As in the case of a 

single qubit QND interaction with bath 0, Hij|, the density matrix is evaluated in the system eigenbasis \i n ) = | ± |) 
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(the possible eigenstatcs of J" with eigenvalues j n = ±|). The system-plus-reservoir composite is closed and hence 
obeys a unitary evolution given, in the interaction picture, by 



p(t) = U I (t)p(0)U}(t), 

where 

U I (t)=Te- m) fo dt ' Hl(t '\ 
with Hj(t) = e^ Hs+HR ^/ h H S Rer^ Hs+Hll ^/ h , and T denotes time ordering. Also 

p(0)=p s (0)p R (0), 

i.e., we assume separable initial conditions. Here 

P s (0) = pl uh (0) ® pt 2>j2 (0) • • • ® (0), 



(6) 

(7) 
(8) 
(9) 



is the initial state of the qubit system with p" m j m (Q) = (im.\p s (0 )\in t)- In Eq. (8J, p_r(0) is the initial density matrix 
of the reservoir which we take to be a squeezed thermal bath [J, [3J, [35[ given by 



PR 



(0) = S(r,$)p th &{r,*), 



where 



fe 



e -/3Sajfe"| — j3hu) k b' k b k 



(10) 



(11) 



is the density matrix of the thermal bath at temperature T, with /3 = l/ffegT), fc^ being the Boltzmann constant, 
and 



S(r k) $ k ) = exp 



/b 2 tf 2 



(12) 



is the squeezing operator with r k , being the squeezing parameters |36| . 

In order to obtain the reduced dynamics of the system , we trace over the reservoir variables. The matrix elements 
of the reduced density matrix in the system eigenbasis are obtained for the independent and collective decoherence 
models as: 



A. Independent decoherence model 



/>{w„}W = cx p[*{©{L (*) - A {i„,j„}(*)}] e M-rt indn}{s ^(t)]pt in , jn} (o). (is) 

Here P| injJ - n }(*) stands for (ij,, ii|Trflp a (i)|jL, Ji-i, • and the symbol {i„,j n } stands collectively for 

iii ji> *2, i2i ii, Jl- The superscript m is to indicate that these expressions are for the independent decoherence 
model and the subscript sq indicates that the bath starts in a squeezed thermal initial state. As seen from the 
expressions given below, 0V" ^ j and A 1 ^ ^ j are independent of the bath initial conditions and are given in the 
continuum limit (assuming a quasi-continuous bath spectrum) by 

OO £ 

6 {L,Jn}(*) = 2 [ duI(u)S(w,t) J2 (Wn-jmjn)cOs(wi s ), (14) 

£ m=l,„=2 



A$ ntjn} (t) = 2j du,I(u;)C(u;,t) Y, 



i m j n sm(u>t s ). 



(15) 
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In the above equations, I(u>) is the bath spectral density which for the Ohmic case considered here has the form 



J(w) = — ue^^% 



(16) 



where 70 and u c are two bath parameters. Also 



and 



C(uj,t) 



Lot — sin(wi) 



1 — cos(wi) 



Also in Eqs. (JHJ) and (JUJ), ui s = fc.r m „ ja3|. In Eq. (I3J) r l "(i) is given as 



(17) 



(18) 



{in J„}(sq) 



(t) 



awl (ui) coth( - ) 



3 h(2r)C(^,i){ £( 



m— 1 



1j 

2 £ (i m - jm)(in - J„)cOs(wt s )| 



1=1,71 = 2 



in 2 ( y ) sinh(2r)| cos(w(i - 2a)) (i m - j m ) 2 cos(wi corr i) 



2 . o.,u)t 
—5 sin 
or 



2 £ (i m - Jm)(»n - in) cos(wt corr2 )] + sin(u(t - 2a)) (i m - j m ) 2 sin(w* corr i) 



TT1=1,T1=2 



1^ 

2 X] (* m ~ J"*)(*n - 3n) Sm(l^ corr2 )] } 



1=l,Tl = 2 



(19) 



where we have dehned two new time scales iot corr \ = 2k.r m and u:t corr 2 = k.(r n + r m ) which are due to the non- 
stationary effects introduced by the squeezed thermal bath. Here wc have for simplicity taken the squeezed bath 
parameters as 



cosh (2r(a;)) = cosh(2r), sinh (2r(w)) = sinh(2r), 
&(u>) = aoj, 



(20) 



where a is a constant depending upon the squeezed bath. 



B. Collective decoherence model 



The reduced density matrix is given by 



col 

{indn} 



m(*)Kw„}(°)- 



(21) 



The superscript col is to indicate that these expressions are for the collective decoherence model and the subscript sq 
indicates that the bath starts in a squeezed thermal initial state. As in the case of independent decoherence, co1 and 
A co1 are independent of the bath initial conditions and are given in the continuum limit (assuming a quasi-continuous 
bath spectrum) by 



00 

\t) = J duI{u)S{w,t) 



(£z m ) 2 -(£j m ) 5 



rn — 1 



m=l 



(22) 
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A col (t) = 0. (23) 



The bath spectral density I(oj) is as in Eq. HI]). In Eq. (H]), ^\f n Jn }( sq )(i) is 

th(^) 

m— 1 
L 



cosh(2r)C(w,t)[^ (i m - ]m )] 2 

l m=l 

— sin 2 (y ) sinh(2r)| cos(w(i - 2a))[^ (i m - j m ) 2 cos(u)t corr i) 



m—1 



2 ^ {im ~ jm)(in ~ jn) cos(ujt corr2 )] + sm(uj(t - 2a))[^(i m - j m ) 2 sin(a;t corr i) 



th = 1,ti = 2 



2 X] (* m ~ 3™){in ~ jn)sin(W corr2 )]} 



1 = 1,71 = 2 



(24) 



All the other terms arc as defined above. On comparing Eq. (|24|) with (|19|) . we find that the terms proportional to 
sinh(2r), arising from the non-stationarity of the squeezed bath, are same while the terms proportional to cosh(2r) 
differ from each other. For the collective decoherence model, ojt s = k.r mn = 0, but the two time-scales coming from 
the non-stationary components of the squeezed thermal bath, i.e., ujt corr \ = 2k.r m and Lut corr 2 = k.(r n + r m ) are 
both non-zero, indicative of correlations induced between the qubits by the bath squeezing. For the case of zero bath 
squeezing, both the Eqs. (j2"4")) and (fT9|) reduce to their corresponding values for the case of a thermal bath [33j]. 

IV. TWO QUBIT INTERACTION 

Here we specialize the general considerations of the previous section to the case of two qubits. 

A. Independent decoherence model 

The reduced density matrix is a specialization of Eq. (fTB")) to the case of two qubits, say a and b. Here p s ^ j , (t) 
would be p% j b }(t) w hich represents (i a: ib\T^ rP s {t)\j ai jb) ! where the states |i a ) or \ib) have eigenvalues ±i. We will 
collectively represent the two-particle index ab by a single 4-level index according to the following scheme: 

11 11 11 11 

— , — =0, --,- = 1, -, — = 2, -,- = 3. 
2' 2 ' 2'2 ' 2' 2 ' 2'2 

Thus there will be sixteen elements of the density matrix, which we enumerate below. They are seen to satisfy the 
symmetries: 

0Sa(f)=P25(*) = Poi(*) = Pio(*), (25) 

where * in the superscript indicates complex conjugation, and of course the first and last equality follow from the 
hermiticity of the density operator. In the Eqs. ([23]) . 9V" - \{t), AV" ■ ,(t) can be obtained from the Eqs. (|14[) . 

(fl"5|) . respectively, and j„}( sq )(*) from the Eq. (fT9|) and are given by 

oo 

Of 2 (t) = Q^t) = [ duI(cj)S{u>,t)cos(ut s ), (26) 



A™(t) = AjJi(t) = - / du;I(u})C(u},t)sm{ujt s ), (27) 
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e&(t) = ©\oW = -e&(*) = -eft(t), 
A&(f) = Ai n (t) = -A*(t) = -AS(t), 



(28) 
(29) 



and 



oo 

™(t) = f dujl{uj) coth(^^) cosh(2r)C*(w, t) - \ sin 2 A sinh(2r){ cos(o;(i - 2a)) cos{ujt [ ^ rrl ) 





sin 



in(u;(f - 2a)) sin(a42. rl )} 



(30) 



for all the above combinations. In the above equations, u>t s stands for k ■ r a t while uit^J rrl = 2k ■ r&. Interestingly, for 
the above cases, the correlation time ojt corr 2 = k.(r a + r&) is absent. It can be seen that 

Plait) = plM, (a -0,1, 2, 3), (31) 

from which follows that the population remains unchanged. This is a consequence of QND nature of the S — R 
interaction. Also, 



P5i(*) = Pi5(*)=Pia(*), 

P03(*)=P03(*)> 

i.e., these components are purely real. In the Eqs. ([32]) . 6 m (i), A ln (i) and r™(i) are given by 

6 in (t) = = A m (t), 



(32) 
(33) 



and 



sq,30 



(t) 



sq,03 



(*) 



/3hu> 

duil{ui) coth(— — ) 



2 o 

2cosh(2r)C(w,t)[l + cos(wt s )l r sin 2 ( — ) sinh(2r) 

w z 2 



x | cos(cj(t — 2a))[cos(2/c • r a ) + cos(2fc • r^) + 2cos(fc • [r a + r&])] 
+ sin(u;(t - 2a))[sin(2fc • r a ) + sin(2fc • r b ) + 2sin(fc • [r a + r 6 ])]| 



(34) 



- sq,21 



(t) 



sq,12 



(*) 



(hoi (lu) coth( - ) 



2cosh(2r)C(u>,t)[l - cos(wi s )] - sin 2 (y ) sinh(2r 



x | cos(u>(t — 2a))[cos(2/c • r a ) + cos(2fc ■ ffe) — 2 cos(fc ■ [r a + rb})} 

+ sin(w(i-2a))[sin(2fc-r a ) +sin(2fc-r 6 ) -2sin(fc- [r a + r b ])]| . (35) 
Thus we see that the Eqs. ([34|). (|35l) . depend on both 2k-r a and 2fc-n,, and iot corr 2 which is as defined above. Further, 

P3lW=^W = P02(0=P20W, (36) 

where * in the superscript indicates Hermitian conjugation. In the Eqs. (f3"u| . ln (t), A la (t) are 



©alW = Soa(*) = / dw/(w)S(w,t)cos(wt s ), 



(37) 



A£(f) = Aj, n 2 (<) = / tiw/(a;)C(w,t)siii(a;t a ), 



(38) 
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FIG. 1: Purity as a function of temperature T (in units where % = Ub = 1) for the independent decoherence model. The 
bold, large-dashed and small-dashed curves correspond to evolution time t = 2.0, 3.0 and 2.0, respectively. The bath squeezing 
parameter r (1201) is equal to 0.2, 0.2 and 1.0, respectively, for the three curves. Here and in all the subsequent figures, the 
squeezing parameter a (|20|) is set equal to zero. Also the bath parameters 70, u) c (|16|) are equal to 0.01 and 100.0, respectively. 
All the inter-qubit distances are defined on the scale of the resonant wavelength coming from the wavevector k ((5} as a result 
of the position dependent couplings of the qubits with the bath. In all the plots concerned with the independent decoherence 
model, kr 12 is set equal to 1.1. 



&m = = -e&(t) = -eis(t), 



(39) 



A%(f) = AjS(t) - -A™(t) = -A£(t), 



(40) 



and 



rjS(t) 



dul(u>) coth(^) 





sin 



cosh(2r)C(o>, t) - ^ sin 2 (y ) sinh(2r){ cos(tj(i - 2a)) cos^i^^) 



(w(t-2a))sm(a;tf ) T . rl )} 



(41) 



for all the above combinations. In the above equations, ojt s stands for k ■ r a t, while ut corrl = 2k ■ r a . Interestingly, for 
the above cases, the correlation time ujt corr 2 = k.(r a + r&) is absent. The Eqs. (|2"5|) . (|5T|) . ([52^1 and (|5r]|) cover all the 
density matrices for the two-qubit independent decoherence model. 

Figure |T|) depicts the behavior of purity, defined here as Tr(p 2 (t)) for p(t) as obtained in this subsection for the 
independent decoherence model, as a function of temperature T for various evolution times and bath squeezing r 
(|20[) . In all the figures in this article, we consider the initial state to be an equal superposition state, which can be 
obtained by applying H ® H on the state |0) = | — 5, —5), where H is the Hadamard transformation. It can be seen 
that with the increase in temperature, as also evolution time t and bath squeezing r, the system becomes more mixed 
and hence looses its purity. 



B. Collective decoherence model 



The reduced density matrix is a specialization of Eq. (|21j) to the case of two qubits, say a and b. The notations 
are as before. 



= pS(t) = p5i(*) = Pio(*). 



(42) 



where * in the superscript indicates Hcrmitian conjugation. In the Eqs. (|4"2")h col (i), A col (i) (= 0) are obtained from 
the Eqs. (HH), ([23]), respectively and T c s °\t) from Eq. j24]). They are given by 



e&(t) = e 



01 



00 

(t) = J dioI{u)S(cu,t), 



(43) 
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©23 (*) = ©TO = -©32 (*) = -©Si (*), (44) 
and r™ (t) is as in Eq. (|3T))) for all the cases in Eq. P2"| . with u)t corr \ and uit corr 2 as defined there. As before. 

Plait) = pL(0), (a -0,1, 2, 3). (45) 
This is indicative of QND nature of the S — R interaction which preserves the population. Also, 

til® = pim = pi 2 (t), 

Pso(*)= =Pos(*)- (46) 
In the Eqs. 9 col (i), A col (i) and r^(t) are given by 

6 col (t) =0 = A col (i), (47) 



and 



^sq,3o(*) — -^sq,03(*) 



duil (u>) coth( - ) 



2 . 2 U)t^ 



4 cosh(2r)CV, t) ~ sin^( — ) sinh(2r) 

w J 2 



o 

X 



I cos(w(i — 2a))[cos(2fc • r a ) + cos(2fc ■ r&) + 2 cos(/c • [r a + r&])] 
+ sin(w(* - 2a))[sin(2fc ■ r a ) + sin(2/c • r b ) + 2sin(/c • [r a + r b ])]\ , (48) 

oo 

= -2 / — 7(cj)coth(-^-)sin 2 (y)sinh(2r) 
o 

x [cos(aj(t — 2a))[cos(2fc ■ r a ) + cos(2fc • rb) — 2 cos(fc ■ [r a + r&])] 

+ sin(w(* - 2a))[sin(2fc • r a ) + sin(2fc ■ r&) - 2 sin(fc • [r a + r b })}} . (49) 

It is interesting to note from Eqs. P6"| . (|T7|) and (j4"9")l . that for the case of a purely thermal bath with zero bath 
squeezing, r™ 1 21 (i) = rg°' 12 (t) = 0, thereby implying that for these cases, the corresponding density matrix elements 
do not decay even though they are interacting with the bath. Also, since in a QND S — R interaction, the diagonal 
terms p\ t \ and p 2j2 do not change, this implies that any state a\l) + /3|2) in the subspace span {|1),|2)} remains 
invariant, thereby leading to a decoherence-free subspace. This behavior is not retained, in general, for the case of a 
squeezed thermal bath. However, since a thermal bath is the asymptotic limit of a squeezed thermal bath [37l. l38l. [39| . 
decoherence-free subspaces emerge asymptotically. Further, 

Ph(t) = Pis® = Pot® = P*2o(t), (50) 
where * in the superscript indicates complex conjugation. In Eq. ([50]) , 8 col (£), A col (i) (= 0) are 



©si (*) = ©S°2 (*) = J duI(fj)S(u,t), (51) 



Qtf(t) = 02oW = -e§?(t) = -e c f(t), (52) 

and T™'(i) is as in Eq. (001), with 2k ■ r b —> 2k ■ r a . The Eqs. l[4"2"l). (|i5]l. (gU) and dSUJ) cover all the density matrices 
for the two-qubit collective decoherence model. 

Figure ([2]) depicts the behavior of purity, defined here as Tr(p 2 (i)) for pit) as obtained in this subsection for the 
collective decoherence model, as a function of temperature T for various evolution times and bath squeezing r (|20l) . 
As before, it can be seen that with the increase in temperature, as also evolution time t and bath squeezing r, the 
system becomes more mixed and hence looses its purity. 
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FIG. 2: Purity as a function of temperature T (in units where fi = ks — 1) for the collective decoherence model. The bold, 
large-dashed and small-dashed curves correspond to evolution time t = 2.0, 3.0 and 2.0, respectively. The bath squeezing 
parameter r (|20|) is equal to 0.2, 0.2 and 1.0, respectively, for the three curves. All the inter-qubit distances are defined on the 
scale of the resonant wavelength coming from the wavevector k ([5| as a result of the position dependent couplings of the qubits 
with the bath. In all the plots concerned with the collective decoherence model, kr\i is set equal to 0.05. 



V. SYMMETRIES IN THE DYNAMICAL SYSTEM 



In this section we consider the twoqubit evolution, developed in the previous section from the point of view of some 
fundamental symmetries. This enables us to view the dynamics from a fresh perspective and is also interesting from 
its own point of view. 

Employing the two-particle index notation used in the previous section, we find that the transformation connecting 
the initial and final density operations can be given by the following operation 



p a b(t) = C ab (t)p ab (0). 



(53) 



The non-trivial aspect of the dynamics that this relation represents is that here C represents, not a matrix, but a 
two-dimensional array, and the multiplication is done element-wise. 

The most general array C that satisfies the symmetry, following only from the hermiticity of p(0) and p(t) is: 



L = 



l a c 2 c 3 

c* 1 c 4 c 5 

c 2 C 4 1 C 6 



(54) 



However, further constraints on the structure of C appear because the dynamical evolution due to QND interaction 
respects spin- flip symmetry (see Eq. (|59p below), which is for example (given for clarity, in the single-qubit notation): 



i.i i = C_i i. 

2 ' 2 ' 2 2 ' 2 ' 



This has the effect that C3 and C4 are real, which we denote by n and r 2 respectively. Further c 2 = 
These are seen by noting that: 



(55) 

5 and ci = eg. 



C3 


= r 

— i -e,g 




— r 


- >~e,g 


= c 3> 


Ci 


= r 


'-'a,s 


— r 


— C* 


= 


<"2 


= r 




— r 

— J -'g,s 




= c 5. 


Cl 


= r 




— r 


— c* 


= c 6, 



(56) 

where the first and third equalities in each equation follow from hermiticity. Accordingly, Eq. (|54|) can be rewritten 



C 



/ 1 ci c 2 r\ 

c* 1 r 2 c* 2 

c* 2 r 2 1 c\ 

V n c 2 ci 1 



(57) 
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(a) (b) 

FIG. 3: Concurrence C (|60[) as a function of time of evolution t at T = 5.0 and bath squeezing parameter r (|20[1 equal to 0.2. 
Figure (a) refers to the independent decoherence model and (b) the collective decoherence model. 



Consider the operator C corresponding to Cj t k, defined by: 

C = J2^.k\j)(kl (J, k = 0,1,2,3). (58) 
Oik 

The spin-flip symmetry can be represented by 

Y, = <j x ®a x . (59) 

Since X = (— icr x ) ® (i<J x ), the above spin- flip symmetry may be described as a rotational symmetry, with angle ir/2 
(resp. 3n/2) in the first (resp. second) qubit coordinate about the a:-axis. 



VI. ENTANGLEMENT ANALYSIS 



In this section, we will study the development of entanglement in the two qubit system, both for the independent 
as well as the collective decoherence model. A well known measure of MSE is the concurrence [27| defined as 

C = max(0, VAi- \/% - \/% - \/%), (60) 
where Xi are the eigenvalues of the matrix 

R = pp, (61) 

with p = (jy ® a y p*a y ® a y and a y is the usual Pauli matrix. C is zero for unentangled states and one for maximally 
entangled states. 

In figure ([3] (a)), we plot the concurrence (|60|) with respect to time for the case of the independent decoherence 
model, while figure ([3] (b)) depicts the temporal behavior of concurrence for the collective decoherence model. It 
is clearly seen from the figures that the two qubit system is initially unentangled, but with time there is a build 
up of entanglement between them as a result of their interaction with the bath. Also the entanglement builds up 
more quickly in the collective decoherence model when compared to the independent model. This is expected as the 
effective interaction between the two qubits is stronger in the collective case. 

Now we take up the issue of entanglement from the perspective of the PDF as in Eq. ((3|). In figures (jH (a)) and 
(b), we plot the weights u>\, u>2, UJ3 and lo^ ^ of the entanglement densities of the projection operators of the various 
subspaces which span the two qubit Hilbcrt space with respect to T for the independent and collective decoherence 
models, respectively. Since lo\ is the weight for the one dimensional projection and hence a pure state and W4 that 
of the maximally mixed state with loi and U3 being intermediary, these plots depict the variation in the contribution 
of the various subspaces to the entanglement of the two qubit system as T increases. As can be seen from both the 
figures, with increase in temperature T, the weight uj\, depicting the pure state component monotonically decreases, 
while the other weights start from zero at T — and increase. Eventually, the weight lo^ depicting a maximally mixed 
state would be expected to dominate, though for the parameter range used in the plots, this feature is not seen. This 
feature of the dynamics of the reduced two-qubit system, specially in the case of the collective decoherence model, 
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FIG. 4: The weights J3| as a function of T, with an evolution time t = 5 and bath squeezing parameter r (|20[) equal to 0.2. 
Figure (a) refers to the independent decoherence model and (b) the collective decoherence model. In both the figures, the bold 
curve corresponds to the weight cji, while the large-dashed, small-dashed and dotted curves correspond to the weights U2, 
and u>4, respectively. 

P 4 (S) 




0.2 0.4 0.6 0.8 1 



FIG. 5: The PDF for the four dimensional projection spanning the full Hilbert space THJI4). 



has an interesting application which will be discussed in detail in Section VII where it will be seen to obey an effective 
temperature dependent Hamiltonian, bringing out the persistence of entanglement even at finite temperatures. In the 
case of the collective decoherence model, the weights loi and LU3 have a greater growth than that for the independent 
decoherence model, depicting the greater entanglement development in the collective model as is also borne out by 
the concurrence plots in figure 

As explained in Section II, the characterization of MSE for a two qubit system via the PDF involves the distri- 
bution functions of four projection operators, III, II2, II3, II4 corresponding to one, two, three and four dimensional 
projections, respectively. These will be represented here as Vi(£), Vi{£\ ^M^-) an d Vi{£), respectively Also, as 
discussed above, Vi(£) would be universal for the two qubit density matrices and would involve the Haar measure on 
5(7(4) j4fj. This is depicted in figure ((5]) and is common to all the two qubit PDF of entanglement. 

Now we consider the T , 2{£) and 7^3 (£) density functions for some representative states of the two qubit system, 
both for the independent as well as collective decoherence models. This enables us to compare the entanglement in 
the respective subspaces of the system Hilbert space. We also plot the full entanglement density function curve V(£) 
with respect to the entanglement £ , at a particular time t. This will enable us to look at the contribution to the 
entanglement from the different projections. 

Figures ([5] (a)) and (b) depict the behavior of the density function V?,{£) for the bath evolution time t = 2.0 
and T = 5.0 for the independent and collective decoherence models, respectively. For these conditions, the value 
of concurrence (|60[) is 0.18 for the case of the independent decoherence model and 0.27 for the collective model, 
depicting the greater entanglement content in the later compared to the former. This is also borne out by these 
figures. As shown in [30j ] . the concurrence for a two dimensional projection is Cn 2 = {£max ~ £cusp)/2. Thus while 
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the figure ([6] (a)) has £ m ax = £cusp, resulting in Cn 2 =0, it is finite for the case of figure ([6] (b)) , thereby reflecting 
the greater entanglement content in the collective model. Figures ([7J (a)) and (b) depict the behavior of the density 
function V^S) for the bath evolution time t = 2.0 and T = 5.0 for the independent and collective decoherence models, 
respectively. In figure ([7] (a)), the parameter £± is w 0.19 while it is one for the collective model as depicted in figure 
(UJ (b)). Since £j_ is evaluated from the state perpendicular to the third non-separable basis in the canonical basis of 
a three dimensional projection [3Cj |. this clearly brings out the greater entanglement content in the collective model 
as compared to the independent case. In figures ([5] (a)) and (b), we depict the full density function V{£ ) for the 
independent and collective decoherence models, respectively and with the same parameters as above. This is obtained 
by a weighted sum over all the contributions of the projection operators from the different subspaces ([3]). Here and in 
all the subsequent figures for the full PDF, the contribution from the one dimensional projection n 1; which is a delta 
function, is represented by a line of height equal to its weight ^ and the point on the absicca is determined by its 
corresponding entanglement. The figures clearly indicate, that for the given bath parameters and evolution time, the 
system is still possesed with considerable entanglement, which is clearly greater in the collective case because of the 
contribution to the full PDF coming almost exclusively from the entanglement in the one dimensional projection. 

Figures (0 (a)) and (b) exhibit the behavior of the density function Vi{£) for the bath evolution time t = 10.0 
and T = 5.0 for the independent and collective decoherence models, respectively This thus depicts a situation where 
the two qubit system has been exposed to the bath for a longer time. For these conditions, concurrence (|60[) is 0.17 
for the case of the independent decoherence model and 0.35 for the collective model, showing greater entanglement 
content in the later compared to the former and also indicating a builtup of entanglement with time. As in the 
corresponding case before, figure ([S] (a)) has £ max = £ C usp, resulting in Cn 2 = 0, while it is finite for the case of figure 
(jni(b)) , thereby reflecting greater entanglement content in the collective model. Figures (fTUl (a)) and (b) depict the 
behavior of the density function V3(£) for the bath evolution time t = 10.0 and T = 5.0 for the independent and 
collective decoherence models, respectively. In figure (fTUl (a)), the parameter £± is w 0.79 while it is one for the 
collective model as in figure (|10l (b)) bringing out greater entanglement content in the collective model in comparison 
with the independent case. In figures (fTTI (a)) and (b), the full density function V{£ ) is given for the independent 
and collective decoherence models, respectively and with the same parameters as above. The point on the absicca, 
for the one dimensional projection, can be clearly seen to have shifted to the right as compared to the corresponding 
case in the figures ([8] (a)) and (b), thereby indicating a builtup of entanglement with time. Also when compared to 
figures ((5]( a )) and (b), the full PDF as in figures QUI (a)) and (b) exhibit a richer structure which is a consequence of 
the additional contribution to the full PDF from the two and three dimensional projections. The fourth dimensional 
projection has an insignificant contribution for the given parameters. 

Figures (fT2l (a)) and (b) give the full density function V{£) for the independent and collective decoherence models, 
respectively with a bath evolution time t = 10.0 and T = 50.0. For these conditions, the value of concurrence ([50]) 
is 0, which would indicate a complete breakdown of entanglement. This would be expected as with the increase in 
the bath temperature T, the effect of entanglement would be destroyed quickly. This is partially borne out by the 
fact that for this case Cn 2 = 0. However, as seen from figure (fT2l (b)), the PDF for the full density function still 
exhibits a rich entanglement structure, coming principally from the contributions from the one and three dimensional 
projections. In contrast , figure (|12l fa)), for the independent decoherence model, exhibits the Harr measure on SU(4) 
and thus represents a maximally mixed state. 

Figure (|13p represents the full density function V(£) for the independent decoherence model with an evolution 
time t = 10.0, T = 20.0 and bath squeezing parameter r equal to 0.2. This would be analogous to that discussed 
in [30l | for NMR quantum computation where concurrence would be zero, and the excess of entangled states over 
the unpolarized background (exhibited by the uniform distribution coming from the density function related to 
the four dimensional projection) is exploited as a resource allowing for non-trivial gate operations, thus depicting 
pseudopure states over the four dimensional background, with the excess being the "deviation density matrix" . 

VII. EFFECTIVE TEMPERATURE DEPENDENT DYNAMICS IN THE COLLECTIVE 
DECOHERENCE MODEL: A BRIEF DISCUSSION 

In a QND S — R interaction, the reduced density matrix of the system does not approach a unique distribution 
asymptotically [J. As seen from figure (fT2l (b)), the PDF for the full density function (for the collective decoherence 
model) exhibits a rich entanglement structure, coming principally from the contributions from the one and three 
dimensional projections which carry equal weights. This feature is seen to persist for higher temperatures and evolution 
times, for the collective decoherence model, with the weights of the subspaces spanned by the four projection operators 
of the PDF remaining intact. From this emerges the fact that for the collective decoherence model, studied here, as the 
effect of the bath on the system increases, the PDF instead of becoming uniform, as expected, gets distributed between 
the subspaces spanned by the one and three dimensional projection operators suggesting a tendency of the system to 
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FIG. 6: The density function V<2,{£} with respect to the entanglement £ for an evolution time t = 2.0, T — 5.0 and bath squeezing 
parameter r equal to 0.2. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence model. 
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FIG. 7: The density function Vs{£ ) with respect to the entanglement £ for an evolution time t = 2.0, T = 5.0 and bath squeezing 
parameter r equal to 0.2. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence model. 

resist randomization. Such a state of affairs would be encountered if the effect of the bath is not a counterpart of the 
collision term (in a Boltzmann equation) , but is more like a Vlasov term, causing long range mean field contributions 

a 

From the numerical results, it is seen that the effect of the bath can be mapped to a T dependent effective 
hamiltonian whose energy eigenvalues scale with temperature. The cigenstates are given by the standard Bell states 
with the ground state being \B\) = 4j(|01) + 1 10) ) while the orthogonal singlet state (|B 4 j = ^(|01) - 1 10) )) is the 
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FIG. 8: The full density function V(£) (|3| with respect to the entanglement £ for an evolution time t = 2.0, T = 5.0 and 
bath squeezing parameter r equal to 0.2. Figure (a) refers to the independent decoherence model and (b) to the collective 
decoherence model. 
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FIG. 9: The density function T- , 2{£) with respect to the entanglement £ for an evolution time t — 10.0, T = 5.0 and bath 
squeezing parameter r equal to 0.2. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence 
model. 
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FIG. 10: The density function Vz(£) with respect to the entanglement £ for an evolution time t = 10.0, T = 5.0 and bath 
squeezing parameter r equal to 0.2. Figure (a) refers to the independent decoherence model and (b) to the collective decoherence 
model. 



highest energy state, and is practically decoupled (with no population). The next excited state is degenerate, with 
two Bell states QB2) = ( |00) — |11)), I-B3) = -1=(|00) + 1 1 1) ) ) spanning the two dimensional subspace. The effective 
temperature dependent hamiltonian is given by 

4 

H cfi = Y / E l (P)\B t )(B l \, 

i=l 

where \B X ), \B 2 )(B 2 \ + \B 3 )(B 3 \ = |00)(00| + |H)(H| and \B 4 ) are the Bell states, as defined above, with eigenvalues 
Ai = 0.5, A2 = A3 = 0.25 and A4 as 0, respectively. Since the Bell states are completely entangled, the effective 
hamiltonian has no linear terms in the qubit polarizations and has the form 

Heff ~ ^(1 - (-*iM 2) + - <W)) + + 2^^), (62) 

in writing which the singlet term has been dropped, as it is energetically very far separated from the other three levels. 
The above analysis places in perspective the surprising result that although the system is evolving, through the effective 
Hamiltonian, the entanglement density function remains practically restricted to the 3-dimcnsional subspace, with 
a large contribution from a Bell state, as a signal to the 3-dimensional background. The restriction of the effective 
dynamics from four to three levels is also seen in the case of two qubit evolution via a dissipative S — R interaction 
with a thermal bath initially at T = j^, for the collective decoherence model. However, there the reason for it is 
simply given by the fact that for the above conditions, the coupling term connecting one of the levels to the others 
goes to zero, thereby reducing the dynamics to that between three levels. 

An interesting analog of the discussion in this Section comes in the work presented in [43[ . There it was shown by 
the authors that for the scenario where there exists a system consisting of three subsystems with the first and the third 
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FIG. 12: The full density function V(E) © with respect to the entanglement S for an evolution time t — 10.0, T = 50.0 and 
bath squeezing parameter r equal to 0.2. Figure (a) refers to the independent decoherence model and (b) to the collective 
decoherence model. 



interacting with each other via their interactions with the mediating second subsystem, a signature of entanglement 
between the first and the third subsystems is the degeneracy in the ground state of the system. Here we have a similar 
situation with the two qubits interacting with a bath which in turn mediates the inter-qubit interaction. From our 
effective Hamiltonian H e ff, we see that the first excited state (not the ground state), spanned by the Bell states \B 2 ) 
and I-B3), is degenerate and the system exhibits a strong entanglement even at finite temperatures. Another work 
by the same authors |44j |. studied the persistence of mixed state entanglement at finite T. This would be important 
as quantum effects can be expected to dominate in regions where entanglement is nonzero. They considered the 
transverse Ising model and studied the two-site entanglement, using concurrence as the entanglement measure, and 
found appreciable entanglement in the system at finite T above the ground state energy gap, one of their motivations 
being the influence of nearby critical points to the finite T entanglement. The persistence of entanglement in a two- 
qubit system interacting with the bath via a purely dephasing interaction (QND) would suggest a broad applicability 
of these concepts, thereby highlighting the interconnection of ideas of quantum information to quantum statistical 
mechanics. 



VIII. AN APPLICATION TO QUANTUM COMMUNICATION: QUANTUM REPEATERS 

We now make an application of the two-qubit reduced dynamics obtained from QND system-reservoir interaction to 
a quantum repeater [3 if . used for quantum communication over long distances. The efficiency of quantum communi- 
cation over long distances is reduced due to the effect of noise, which can be considered as a natural open system effect. 
For distances much longer than the coherence length of a noisy quantum channel, the fidelity of transmission is usually 
so low that standard purification methods are not applicable. In a quantum repeater set-up, the channel is divided 
into shorter segments that are purified separately and then connected by the method of entanglement swapping, which 



17 




is the quantum teleportation [4a ] of entanglement. This method can be much more efficient than schemes based on 
quantum error correction, as it makes explicit use of two-way classical communication. The quantum repeater system 
allows entanglement purification over arbitrary long channels and tolerates errors on the percent level. It requires a 
polynomial overhead in time, and an overhead in local resources that grows only logarithmically with the length of 
the channel. 

Here we consider the effect of noise, introduced by imperfect local operations that constitute the protocols of 
entanglement swapping and purification [32], on such a compound channel, and how it can be kept below a certain 
threshold. The noise process studied is the one obtained from the two-qubit reduced dynamics via a QND interaction, 



instead of the depolarizing noise considered in [3l|. A detailed study of the effect of the two-qubit noise on the 
performance of a quantum repeater is underway and will be reported elsewhere. Here we treat this problem in a 
simplified fashion, and study the applicability and efficiency of entanglement purification protocols in the situation of 
imperfect local operations. 

A quantum repeater involves the two tasks of entanglement swapping, involving Bell-state measurements, and 
entanglement purification, involving CNOT gates. The Bell-state measurement may be equivalently replaced by a 
CNOT followed by a projective single-qubit measurement. In entanglement swapping, two distant parties initially not 
sharing entanglement with each other, but sharing entanglement separately with a third party, become entangled by 
virtue of a multi-partite measurement by the third party on the latter's two halves of entanglement. Entanglement 
purification involves two parties employing local operations and classical communication (LOCC) to improve the 
fidelity F of Einstcin-Podolsky-Rosen (EPR) pairs they share, with respect to a maximally entangled state. The local 
operations involve two-qubit gates such as the CNOT operation, followed by single qubit measurement, and a possible 
discarding of an EPR pair. Provided F > 0.5, and at the cost of losing shared (impure) entanglement, the two parties 
can increase the fidelity of the remaining shared entanglement to 

pi = F 2 + [(1 - ^)/3] 2 (63) 

f 2 + [2F(l-F)/3] + (5/9)(l-F) 2 ' k ' 

where F and F' are, respectively, the input and output fidelities of the entanglement purification protocol proposed 
by Bennett et al. [32| . 

In the simplified scenario we consider, the output of the noisy CNOT is taken to be a mixed separable state, in 
place of a pure separable state that is obtained in the noiseless case. As a further simplification, in order to facilitate 
an easy connection with the purification protocol due to Bennett ct al., this mixed state is assumed to be of the form: 

p(F) = F 2 |+, +)(+, +| + (1 - F 2 )|-, -)<-, -|, (64) 

where we use the notation |±) = (l/\/2)(|^) ± |— i)), and |±,±) = |±) ® |±). Thus, p(F) is a mixture in the 
two dimensional space spanned by {| + ,+),| — ,— )}, parametrized by fidelity F, given by (+, + \p(F)\+, +). In the 
notation of Section WK\ these basis states are ±(|0) + |1) + |2) + |3)) and ±(|0) - |1) - |2) + |3)). The state p(F) is 
then the input to the purification protocol, whereby we obtain the output fidelity F' as a function of F. 
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FIG. 14: Purification loop for connecting and purifying EPR pairs. The noisy channel is modelled as a two-qubit QND 
interaction with the environment in the independent (kr a b > 1) (plot (a)) and collective (kr a b <S 1) (plot (b)) decoherence 
regime, with the input state given by p(F) in Eq. (|64p . where r a t is the inter-qubit distance. The bold line is the F = F' 
plot, the small and large dashed curves represent T — 2 and bath squeezing parameter (Eq. ([20}) r = 0.5, with t = 3 and 5, 
respectively. The dotted line represents t = 5, T = 2, r = 1.5, while the dot-dashed curve is due to the noiseless Bennett et al. 
protocol. 



We depict in figures (|T3T a)) and (b) the modified 'purification loop', obtained by subjecting the noiseless loop to 
the above model of noise. In all cases in both figures, as expected, we find that noise degrades the performance of 
the purification protocol. To evaluate the performance, it may be compared with the noiseless case, given by the 
dot-dashed curve in both figures. That this curve lies above the F = F' line in the closed range [0.5, 1] implies 
that fidelities above the minimum value F m in = 0.5 can be corrected to the maximum value i^max = 1 by repeated 
application of purification. The degrading effect of noise can be seen in two ways: it introduces an off-set, whereby 
an input of F = 1 does not yield the same output and restricts F max to values less than 1. 

In both figures, comparison of the small and large dashed curves shows that increasing bath exposure time degrades 
the fidelity by suppressing F max (the point where a curve cuts the F = F' line from above) and increasing the off-set 
at F — 1. A comparison of large-dashed and dotted curves brings out a similar degrading effect of increasing bath 
squeezing. A similar effect can be shown for temperature. The surprise that emerges in comparing the collective 
and independent decoherence models is that, contrary to expectation, the latter produces a less damaging effect on 
purification than the former. 



IX. CONCLUSIONS 



In this article, we have analyzed in detail the dynamics of entanglement in a two-qubit system interacting with its 
environment via a QND S — R interaction. The system and reservoir are initially assumed to be separable with the 
reservoir being in an initial squeezed thermal state. Since the resulting dynamics becomes mixed, in order to analyze 
the ensuing entanglement, we have made use of a recently introduced measure of mixed state entanglement via a 
PDF. This enables us to give a statistical and geometrical characterization of entanglement. 

After developing the general dynamics of N qubits interacting with their bath (reservoir) via a QND S — R 
interaction, we specialized to the two-qubit case for applications. Due to the position dependent coupling of the 
qubits with the bath, the dynamics could be naturally divided into an independent and collective decoherence regime, 
where in the collective decoherence regime, the qubits are close enough to feel the bath collectively. We analyzed the 
open system dynamics of the two qubits, both for the independent as well as the collective regimes and saw that in 
the collective regime, there emerges the possibility of a decoherence-free subspace for the case of zero bath squeezing. 
Interestingly, the dynamics was found to obey two symmetry operations: hermiticity and spin-flip symmetries. The 
existence of the nontrivial spin-flip symmetry would explain the emergence of a decoherence-free subspace (DFS) 4G:j, 
thereby providing a concrete instance of a DFS. 

We then made an analysis of the two-qubit entanglement for different bath parameters. We analyzed both con- 
currence as well as the PDF by finding the entanglement content of the various subspaces that span the two-qubit 
Hilbcrt space. Although the PDF agrees in its predictions with concurrence, it is able to extract more information 
out of the system as a result of its statistical-geometrical nature. Thus we were able to consider an example analogous 
to NMR quantum computation, wherein the concurrence would be zero, and the excess of entangled states over the 
unpolarized background is exploited as a resource allowing for non-trivial quantum information processing. For the 
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collective decoherence model the PDF for the full density function exhibits a rich entanglement structure, coming 
principally from the contributions from the one and three dimensional projections which carry equal weights thereby 
suggesting a tendency of the system to resist randomization. This feature is seen to persist even for higher temper- 
atures and evolution times with the weights of the subspaces spanned by the four projection operators of the PDF 
remaining intact, thereby enabling us to give an effective T dependent dynamics in the collective decoherence regime. 
A comparison of this with some related works suggests the applicability of quantum information theoretic ideas to 
quantum statistical mechanical systems. Finally we made an application to a simplified model of a quantum repeater, 
which can be adapted for quantum communication over long distances. 

Acknowledgments 

We wish to thank Shanthanu Bhardwaj for numerical help. 



[1] W. H. Louisell, Quantum Statistical Properties of Radiation (John Wiley and Sons, 1973). 
[2] A. O. Caldeira and A. J. Leggett, Physica A 121, 587 (1983). 

[3] W. H. Zurek, Phys. Today 44, 36 (1991); Prog. Theor. Phys. 87, 281 (1993) 

[4] S. Banerjee and R. Ghosh, J. Phys. A: Math. Theo. 40, 13735 (2007); eprint |quant-ph/0703054| 
[5] V. B. Braginsky, Yu. I. Vorontsov and K. S. Thorne, Science 209, 547 (1980). 

[6] V. B. Braginsky and F. Ya. Khalili, in Quantum Measurements, edited by K. S. Thorne (Cambridge University Press, 
Cambridge, 1992). 

[7] D. F. Walls and G. J. Milburn, Quantum Optics (Springer, Berlin, 1994). 

[8] W. H. Zurek, in The Wave-Particle Dualism, edited by S. Diner, D. Fargue, G. Lochak and F. Selleri (D. Reidel Publishing 
Company, Dordrecht, 1984). 

[9] C .M. Caves, K. D. Thorne, R. W. P. Drever, V. D. Sandberg and M. Zimmerman, Rev. Mod. Phys. 52, 341 (1980). 
[10] M. F. Bocko and R. Onofrio, Rev. Mod. Phys. 68, 755 (1996). 

[11] C. J. Myatt, B. E. King, Q. A. Turchette, C. A. Sackett, et al, Nature 403, 269 (2000). 

[12] Q. A. Turchette, C. J. Myatt, B. E. King, C. A. Sackett, et al, Phys. Rev. A 62, 053807 (2000). 

[13] G. J. Pryde, J. L. O'Brien, A. G. White, et al, Phys. Rev. Lett. 92, 190402 (2004); J. L. O'Brien, G. J. Pryde, A. G. 

White, et. al, Nature 426, 264 (2003). 
[14] R. Onofrio and L. Viola, Phys. Rev. A 58, 69 (1998). 
[15] J. S. Bell, Physics 1, 195 (1964). 

[16] M. Nielsen and I. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cambridge, 
2000). 

[17] P. Shor, SIAM Journal of Computing 26, 1484 (1997); L. K. Grover, Phys. Rev. Lett. 79, 325 (1997). 

[18] S. Calderbank and P. Shor, Phys. Rev. A 54, 1098 (1996); A. Steane, Proc. Roy. Soc, London, Ser. A 452, 2551 (1996). 

[19] A. Beige, D. Braun, B. Tregenna and P. L. Knight, Phys. Rev. Lett. 85, 1762 (2000). 

[20] M. Fleischhauer, S. F. Yelin and M. D. Lukin, Opt. Commun. 179, 395 (2000). 

[21] S. Schneider and G. J. Milburn, Phys. Rev. A 65, 042107 (2002). 

[22] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press 2002). 
[23] H. Carmichael, An Open Systems Approach to Quantum Optics (Springer 1993). 

[24] M. P. Almeida, F. de Melo, M. Hor-M eyll, A. Salles, et al, Science 316, 579 (2007); A. Salles, F. Melo, M. P. Almeida, 

M. Hor-Mevll. et al, larXiv:0804.455"6l ' 

[25] R. Srikanth and S. Banerjee, Phys. Rev. A 77, 012318 (2008); larXiv:0707.0059l 

[26] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin and W. K. Wootters, Phys. Rev. A 54, 3824 (1996). 
[27] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[28] F. Mintert, A. R. R. Carvalho, M. Kus and A. Buchleitner, Phys. Reports 415, 207 (2005). 
[29] R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

[30] S. Bhardwaj and V. Ravishankar, Phys. Rev. A 77, 022322 (2008). 

[31] H.-J. Briegel, W. Diir, J. I. Cirac and P. Zoller, Phys. Rev. Lett. 81, 5932 (1998); W. Diir, H.-J. Briegel, J. I. Cirac, and 

P. Zoller, Phys. Rev. A 59, 169 (1999). 
[32] C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher et al, Phys. Rev. Lett. 76,722 (1996). 

[33] J. H. Reina, L. Quiroga and N. F. Johnson, Phys. Rev. A 65, 032326 (2002) 

[34] S. Banerjee, J. Ghosh and R. Ghosh, Phys. Rev. A 75, 062106 (20 07); eprint |quant-ph/0703055| 
[35] S. Banerjee and R. Srikanth, Eur. Phys. J. D 46 335 (2008); eprint |quant-ph/061116"Tj ~ 

[36] C. M. Caves and B. L. Schumacher, Phys. Rev. A 31, 3068 (1985); B. L. Schumacher and C. M. Caves, Phys. Rev. A 31, 
3093 (1985). 

[37] B. L. Hu and A. Matacz, Phys. Rev. D 49, 6612 (1994). 
[38] S. Banerjee, Physica A 337, 67 (2004). 



20 



[39] S. Banerjee and J. Kupsch, J. Phys. A: Math. Gen. 38, 5237 (2005). 

[40] T. E. Tilma, M. Byrd and E. C. G. Sudarshan, J. Phys. A: Math. Gen. 35, 9255 (2002). 

[41] E. Carlen, R. Esposito, J. L. Lebowitz, R. Marra and A. Rokhlenko, Preprint (1994) (unpublished). 

[42] Z. Ficek and R. Tanas, Physics Reports, 372, 369 ( 2002). 

[43] H. L. Haselgrove, M. A. Nielsen a nd T. J. Osborne, |arXiv:quant-ph/0308083| 
[44] T. J. Osborne and M. A. Nielsen, |arXiv:quant-ph/0202162"r ~ 

[45] C. H. Bennett, G. Brassard, C. Crepeau, R. Josza, A. Peres and W. K. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 

[46] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 (1998); L. M. Duan and G. C. Guo, Phys. Rev. A 57, 737 (1998); D. 
A. Lidar, I. L. Chuang and K. B. Whaley, Phys. Rev. Lett. 89, 2594 (1998); D. A. Lidar and K. B. Whaley, Irreversible 
Quantum Dynamics, F. Benatti and R. Floreanini (eds.) Springer Lecture Notes in Physics 622 (Berlin 2003). 



